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We present a search for associated Higgs boson production in the process pp — > W/ZH — > £ £ + 
X in ee, e/i, and /i/i final states. The search is based on data collected by the DO experiment at the 
Fermilab Tevatron Collider at yfs — 1.96 TeV corresponding to 5.3 fb _1 of integrated luminosity. 
We require two isolated leptons (electrons or muons) with the same electric charge and additional 
kinematic requirements. No significant excess above background is observed, and we set 95% C.L. 
observed (expected) upper limits on ratio of the production cross section to the standard model 
prediction of 6.4 (7.3) for a Higgs boson mass of 165 GeV and 13.5 (19.8) for a mass of 115 GeV. 

PACS numbers: 14.80.Bn, 13.85. Qk, 13.85.Rm 



INTRODUCTION 



In the standard model (SM), the Higgs boson decays 
predominantly to a WW pair for Higgs boson masses 
above 135 GeV and, with a moderate branching frac- 
tion, to a tt pair for lower masses, both of which de- 
cay to leptonic final states involving neutrinos. Con- 
sequently the associated production of a Higgs boson, 
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pp -)■ W/ZH -> £ ± l± + X, which has additional fi- 
nal state particles, provides easily-detected experimental 
signatures comprising two leptons of the same electric 
charge. This requirement rejects SM processes with op- 
positely charged dilcptons that occur with high produc- 
tion rates such as Z/j*, WW, and tt. Therefore, the like 
charge signature from associated vector boson-Higgs pro- 
duction has an advantage over direct Higgs production, 
pp — > H — > WW , where only unlike charged leptons are 
produced in the final state. 

The DO Collaboration previously published a search for 
associated Higgs production with the like charge dilepton 
signature based on approximately 400 pb _1 of integrated 
luminosity in 2006 [1]. The most recent result from a 
similar search by the CDF experiment using 4.8 fb _1 of 
integrated luminosity [2] was included in the combination 



of Tevatron searches in the H — > WW decay mode in 
2010 [3]. 

In this Article, we present a search for associated 
Higgs boson production with like charged dilcptons using 
5.3 fb _1 of integrated luminosity collected by the DO de- 
tector during the Tevatron Run II period between 2002 
and 2009. The search combines the three leptonic final 
states with either electrons or muons: ee, [i\i and e/z. 

Background processes for the like charge lepton signa- 
tures are diboson production, pp — > W Z — > £v££ and 
pp — > ZZ —$■ ££££. Non-resonant triple vector boson pro- 
duction (VVV, V = W,Z) and the production oiti+V 
are negligible. There are two types of instrumental back- 
grounds. The first, "charge flip" , originates from the mis- 
reconstruction of the lepton charge. For the same lepton 
flavor channels (ee and /i/i) this background arises mainly 
from the Drell-Yan process, pp — > Z/-y* — > £ + £~ . When 
the two leptons are of different flavor, this background is 
negligible. The second instrumental background is from 
falsely identified leptons which originate from jets and 
photons converting to electrons in W boson or multijet 
production. Although these instrumental effects occur 
at low rate, the associated backgrounds make sizeable 
contributions to the dilepton selection due to the large 
production cross sections of the underlying physics pro- 
cesses. 

The Higgs boson signal contains multiple neutrinos in 
the final state, hence a complete reconstruction of the 
Higgs boson mass is not possible. A Higgs boson signal 
would appear as an excess of events with like charged lep- 
tons with kinematic properties consistent with VH decay, 
including missing energy from the neutrinos in the lep- 
tonic decays of the vector bosons, as well as additional fi- 
nal state objects, mainly jets, from the decay of the third 
vector boson. A multivariate technique is employed to 
provide maximum separation between signal and back- 
ground processes based on various kinematic variables. 
In the absence of an excess over the expected number 
of events from background processes, upper cross section 
limits are set. 



EVENT SELECTION 

Events are preselected by identifying at least two lep- 
tons. The selection of the electrons and muons include 
kinematic requirements as well as criteria on their quality. 
The final event selection is performed with a multivariate 
discriminant, which will be discussed later in this Article. 

The DO detector is composed of a central tracking sys- 
tem with a silicon microstrip tracker (SMT) and a central 
fiber tracker (CFT) embedded within a 2 T solenoidal 
magnet, preshower detectors, a uranium/liquid-argon 
calorimeter with electromagnetic (EM) and hadronic sec- 
tions, and a muon spectrometer with drift tubes, scintil- 
lation counters and toroidal magnets [4]. The detector 



was upgraded in spring 2006 to include, among others, an 
additional inner layer of silicon microstrip tracking [5] . A 
lepton is identified by the presence of a track, and its elec- 
trical charge is determined by the direction of the track 
curvature in the magnetic field. Events with like charged 
dilcptons are retained for the analysis while those with 
unlike charged dilcptons are used to validate the event 
reconstruction and simulation. 

Electrons are characterized by their interaction in the 
EM calorimeter and are required to match a track. The 
energy is measured in the EM and the first hadronic 
layers of the calorimeter within a simple cone of radius 
TZ = ^/(Ary) 2 + (A</>) 2 =0.2, where r] and 4> are the pseu- 
dorapidity [6] and the azimuthal angle, respectively. The 
electron cluster must satisfy: calorimeter isolation frac- 
tion, / iso = [E tot {1Z < 0.4)'- E E m(K < 0.2)}/E EM (n < 
0.2), less than 0.2, where E tot is the total energy in 
the isolation cone of radius 1Z = 0.4 and E-em is the 
EM energy in a cone of radius 1Z = 0.2; EM fraction, 
/em = Eftyi/Etot, greater than 0.9, where both energies 
are measured within the cone of TZ = 0.2; and the ra- 
tio of electron cluster energy to track momentum, E /p, 
between 0.5 and 3.0. In addition, the electron candi- 
dates are required to have an eight-variable likelihood 
(£ e ) greater than 0.85, where the likelihood is calculated 
from /i so ; /em; E/p; number of tracks within the isola- 
tion cone; scalar sum of the tracks' transverse momenta, 
Pt, within 0.05 < 1Z < 0.4 of the electron track; track- 
cluster match probability computed from the spatial sep- 
aration and the expected resolution; track distance to the 
primary vertex at closest approach (dca); and, lastly, co- 
variance matrices built from the energy depositions in 
various layers of the calorimeter to represent the longi- 
tudinal and lateral shower development. A loose qual- 
ity electron is defined by relaxing the requirement on the 
likelihood to C e > 0.2, and down to 0.01 or 0.0 depending 
on the purpose, and by removing the E/p requirement to 
estimate the lepton fake rate and model the backgrounds 
from W boson and multijet production. 

Muons are identified by the presence of at least one 
track segment reconstructed in the muon spectrometer 
which is spatially consistent with a track in the central 
detector, where the momentum and charge are measured 
by the curvature of this track. The muon candidate must 
pass cosmic ray veto timing criteria, be outside a cone of 
radius 1Z < 0.1 from any jet of particles present in the 
event, and must not share its track with an electron can- 
didate satisfying the calorimeter isolation and EM frac- 
tion requirements described above. Muon isolation is im- 
posed with two isolation variables defined as the scalar 
sums of the transverse energy, S_Ey al °, in the calorime- 
ter in an annulus with radius 0.1 < 72. < 0.4, and of the 
momenta of tracks, Sp^ k , around the muon candidate 
within a radius of 1Z = 0.5. Each of the two isolation 
variables must be less than 2.5 GeV. Relaxed isolation 
criteria define loose muon quality where the track isola- 



tion upper bound is raised to Sp^ k < 10 GeV and the 
calorimeter isolation is ignored. 

Both leptons are restricted to the central region within 
the CFT coverage, \rj e \ < 1.1 for electrons and |t7 m | < 1.8 
for muons, and have transverse momenta in the range 
15 GeV < pt < 200 GeV, together constraining the 
dilepton invariant mass to 15 GeV < M(£, I) < 250 GeV. 
The upper thresholds were chosen to eliminate imprecise 
charge measurements in the tracking system. The signal 
loss due to the \rj\ requirements is non-negligible, how- 
ever, the reduction of instrumental backgrounds leads to 
an overall gain in the signal-to-background ratio. The 
two lepton tracks are required to have at least one hit 
in the SMT. The longitudinal (transverse) distance sep- 
arating the point of closest approach between the tracks 
and the primary pp interaction vertex must be less than 
1.0 (0.01) cm. The best interaction vertex, which is se- 
lected based on the number of associated tracks and their 
transverse momenta, must be within \z\ < 60 cm, where 
z is the longitudinal coordinate measured from the de- 
tector center. Additional leptons that satisfy the same 
lepton and track quality requirements are allowed up to 
|j7 e | < 2.5 and \t)^\ < 2.0. No upper pt bound is imposed 
for these leptons. 

The hadronic decays of vector bosons and partons 
are identified as jets of particles, and are reconstructed 
by clustering the calorimeter energy deposition within a 
cone of radius 1Z — 0.5 [7]. Jets are required to have 
Pt > 15 GeV and |?j| < 2.5. At least two tracks originat- 
ing from the primary vertex must be found within the 
jet cone. Momentum imbalance in the transverse plane 
of the event implies the presence of the neutrino; the 
missing transverse energy, $t, is reconstructed as the 
negative vector sum of the energies in the calorimeter 
towers and the muon momenta subtracting the calorime- 
ter energy deposition due to muons. 



EVENT SIMULATION 

Associated production of the Higgs boson and the di- 
boson background processes are modeled by Monte Carlo 
(MC) event generators and the detector response is pro- 
vided by the GEANT based [8] simulation of the DO detec- 
tor. The effect of additional pp interactions is reproduced 
by overlaying data events taken from random collisions 
onto the MC generated processes. The signal and dibo- 
son processes are simulated by the PYTHIA event genera- 
tor [9] using CTEQ6L1 [10] parton distribution functions 
(PDFs). The signal event samples are generated for dif- 
ferent Higgs boson masses between 115 GeV and 200 GeV 
with 5 GeV intervals, and normalized to the next-to-next- 
to-leading order (NNLO) cross sections calculated using 
MSTW2008 PDFs [11]. Branching ratios of Higgs boson 
decays are provided by hdecay [12]. The diboson cross 
sections are calculated at next-to-leading order (NLO) 



using CTEQ6.1M PDFs [13], and the transverse momen- 
tum of the diboson system is reweighted to the prediction 
from MC@NLO [14]. Production of W bosons with addi- 
tional emission of partons is simulated with the ALPGEN 
event generator [15], which implements matrix element 
calculations, and the hadronization process is simulated 
by pythia. The NNLO total cross sections calculated 
using MRST2004 PDFs [16] are used to normalize the 
sum of the W production processes with through 5 as- 
sociated light partons and the W production processes 
with heavy flavour quarks (bb or cc). 

The unlike charge dilepton events are used to study the 
event reconstruction and selection as well as instrumental 
effects. The events include physics processes which do 
not contribute to the like charge selection: Z/j*, ti and 
WW production. These processes are simulated using 
the pythia generator, as described above. The Z/-f* 
total cross sections are determined at NNLO, the WW 
cross section at NLO, and the ti cross section is obtained 
at approximate NNLO [17]. 

The overall normalization of the MC samples are ini- 
tially obtained from the corresponding theoretical cross 
sections and the integrated luminosity recorded by the DO 
detector. The kinematic dependences of the efficiencies 
of the triggers, lepton identification and vertex selection 
as well as the mismodcling of instantaneous luminosity 
profile and beam spot position are corrected. Additional 
scale factors specific to this analysis are obtained by com- 
paring the yields of the unlike charge dilepton events 
within the Z resonance in the dilepton invariant mass be- 
tween 70 GeV and 110 GeV. The normalization factors 
include residual effects due to the efficiencies of triggers, 
lepton identification, track quality requirements, vertex 
selection, as well as the measurement of the integrated lu- 
minosity. The efficiencies specific to the electrons and the 
muons are measured with dielectron and dimuon events 
respectively, while the effects that are common to the two 
lepton types have the same multiplicative contributions 
to the total normalization factors in both samples. 



INSTRUMENTAL BACKGROUNDS 



The two largest instrumental backgrounds, charge mis- 
measurement and lepton misidentification, are estimated 
from data, either by measuring their rate using control 
samples enriched with the particular background pro- 
cesses or by performing a fit to the kinematic distribution 
to predict their fraction in the analysis sample. The sep- 
arate contributions to the background from charge flip 
and W and multijet production are discussed below. 



Charge Flip 

The charge flip background, created by mismeasure- 
ment of the charge of one of the leptons, mostly originates 
from the Z/j* process. This occurs when the curvature 
of a high pt track is not correctly measured, or when ad- 
ditional hits from other charged particles and noise are 
present near the track. For electrons, conversion of pho- 
tons from bremsstrahlung radiation is also estimated as 
part of the charge flip background, although the charge 
may be correctly measured. 

Contributions from charge flip in the /x/x sample are 
estimated using two uncorrelated measurements of the 
lepton charge; the first one is the measurement of the 
track curvature in the central tracker, and the second 
measurement, called "local" muon charge, is measured 
by the muon spectrometer. The like charge dimuon sam- 
ples are categorised into three types; events in which the 
two charge measurements give the same answer for both 
leptons (AA), agree for one lepton and disagree for an- 
other one (AD), or disagree for both leptons (DD). The 
number of events in these three categories depends on 
the actual number of like charge [ip events that originate 
from charge mismeasurements, N p , and those that orig- 
inate from true like charge processes such as dibosons, 
jV truo , as well as on the efficiency of the local charge 
measurement, q oc , which gives the probabilities of the 
AA/AD/DD configurations, P, for each case. 
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For a true like charge event to appear in the AA config- 
uration, both leptons must have the correct local charge 
measurements, giving P^A = Qoc ■ tioc- For a charge flip 
event where the central tracker mismeasures the charge 
of one lepton, the local charge must also be incorrect, 
hence P^a = e Zoc(l — tioc)- The parameterization of e; oc 
is measured using unlike charge Z — > /x/z events as a 
function of 1/pt of the muon. The fraction of the charge 
flip events within the like charge /x/x sample, 

/flip = 7V flip /(iV flip + N truc ), 

is determined by solving the over-constrained equations 
yielding / flip = 0.95 ± 0.14 1 . The ratio of the predicted 
number of charge flip events to the number of unlike 
charge events gives a charge mismeasurement rate in the 
dimuon sample of approximately f 0~ 3 . 

A different method is employed to estimate the charge 
flip contribution to the ee background, for which a second 



[1] The positive uncertainty is constrained by the total background 
yield at later stage to obtain the final result. 



charge measurement does not exist. The charge flip event 
rate is measured with data using a control region enriched 
in Z — > e + e~ events, defined by a dielectron invariant 
mass reconstructed with precisely measured energies in 
the calorimeter, 85 GeV < M(e,e) < 100 GeV and an 
azimuthal separation between the leptons of A</>(e, e) > 
2.8. The charge flip rate corresponds to the like charge 
fraction in the Z boson enriched region, removing contri- 
butions from known sources of true like charged dileptons 
such as diboson, W boson and multijet production. The 
measured charge flip event rate in the control region in 
data after the dielectron selection is (8.5 ±1.4) x 10~ 4 . 
The rate is then used to scale the unlike charge distribu- 
tion of the data outside the control region to obtain the 
charge flip prediction for the events selected for analysis. 

The contribution of charge flips to the e/x selection is 
negligible as the dominant Z/"f* production must decay 
via a r lepton pair and is suppressed by the branching 
fraction of the r lepton into an electron or a muon. In ad- 
dition, leptons from t decays have a lower px spectrum, 
hence the average charge flip rate is smaller. 

The kinematic distributions of the charge flip events 
are modeled by unlike charge data, with corrections ap- 
plied for the effects of charge mismeasurement derived 
using MC leptons. The corrections include the depen- 
dence of the charge flip rate on the lepton px and ry in 
the detector, and the resolution of the muon momentum 
measured from the track curvature, parameterized as a 
function of 1/pr- 



W Production 

The production of W bosons contributes to the like 
charge dilcpton background when there is a "fake" lepton 
which can originate from jets or when there are photon 
conversions to electrons. The contribution from these 
backgrounds is estimated using MC simulations, taking 
into account corrections derived from data for the elec- 
tron misidentification rate and photon conversions. 

The corrections to these two subcategories in the MC 
samples are obtained in two steps using template fits to 
the kinematic distributions inaW boson enriched region 
in the ee sample, which requires one of the electrons to 
have loose quality, defined by the L e variable, and to fail 
the tight selection used for the analysis. The first step de- 
termines the normalization of the W boson background 
as well as the charge flip and multijet contributions in 
the control region using a two-dimensional distribution 
defined by M(e,e) and $t- The second step determines 
the fractions of the two contributions within the W bo- 
son sample, jets and photons, while keeping the overall 
normalization of the total sum derived in the first step. 
The presence of a hit in the first layer of the SMT de- 
tector is used to discriminate between the two contribu- 
tions; the photon conversion rate scales as the amount 



of material the photon traverses, hence a sizeable frac- 
tion of electrons from photon conversion are expected 
to have no hit in the first SMT layer. To increase the 
purity of the W process in the control sample, we re- 
quire E T > 25 GeV and M T (e,E T ) > 40 GeV, where 
M T (e,E T ) = y/1-PTe -$t(1-cosA0), where A</> is the 
angle in the azimuthal plane between the electron and 
the direction of $r, and e is the electron which gives the 
maximum value of Mt- The requirement on the electron 
likelihood is relaxed to a minimal value of C e > for the 
first step to maximize the number of W boson events in 
the control region, and raised to 0.01 for the second stage 
to increase the W boson + photon fraction. The normal- 
ization factor obtained in the first stage is 1.2, and the 
correction factors for W + jet and W + photon contri- 
butions obtained in the second step are 0.95 and 1.01, 
respectively. The total systematic uncertainty on the W 
boson background prediction are discussed later. 

The efficiency for a jet identified as a loose electron, 
to pass as a tight electron is determined from a dedi- 
cated dijct data sample where one of the jets is identi- 
fied as a loose quality electron candidate. The events 
are required to contain only one lcpton candidate, back- 
to-back with a jet, A<^)(e,jet) > 2.5, have low transverse 
mass, MT(e,^T) < 40 GeV, and invariant mass with any 
track lower than the Z resonance, M(e,trk) < 60 GeV, 
to suppress contamination from true electrons in W/Z 
decays. Additional physics contamination is estimated 
by comparing the shape of the C e variable for the elec- 
tron candidates in the fake data sample to that in MC 
dijet events. The measured fake rate, parameterized in 
Pt and 77, is then used to scale the W + jet MC sample 
selected with the tight-loose requirement to the tight- 
tight region. No additional correction is made to the W 
+ photon content of the MC sample as the electrons from 
photon conversions are correctly modeled by the simula- 
tion. 



Multijet 

In the case of jets misidcntificd as muons, the multijet 
background contains muons from semileptonic decays of 
heavy flavor quarks, punch-through hadrons in the muon 
detector, and muons from pion or kaon decays in flight. 
In the case of jets misidcntificd as electrons, the multijet 
background contains electrons from semileptonic heavy 
flavor decays, from hadrons misidcntificd as electrons, 
and from photon conversions. 

Multijet contributions to the ee and pp samples are es- 
timated from data events in a control sample containing 
two loose quality like charged leptons using the loose-to- 
tight efficiency of fake leptons, e F , as measured in that 
sample. Within a sample of pure multijet events, the 
fraction of events with no tight lepton (/o) and with ex- 



actly one tight lepton (/i) is 



/o = (l-e F )[l-e F (l-p)}, 
h = 2{l-e F )s F (l-p), 



(2) 



where p is the correlation coefficient, reflecting that the 
identification of a lepton may be different in events where 
another lepton has been already identified. These equa- 
tions are solved for e F and p. Then the fraction of mul- 
tijet events with two leptons passing the tight selection 
criteria is given by 



h 



e F (l - p) + pe F . 



(3) 



The control region is 30 GeV < M(l,l) < 50 GeV 
and A<j>(£,£) > 2.5, where A<j>(£,£) is the azimuthal an- 
gle between the two leptons. The correlation coefficient 
is 0.01 ± 0.05 for the ee channel and 0.18 ± 0.06 for 
the pp channel. The efficiencies are 0.169 ± 0.014 and 
0.091 ±0.009 for the ee and the pp channels, respectively. 
Examination of the variation of e F and p as track qual- 
ity requirements are removed shows no sign of charge flip 
contamination in the control region, and the spectrum of 
Mt(£, $t) shows no sign of W boson contamination. Our 
ability to identify contaminations from other background 
processes (in particular events with charge flips) is lim- 
ited by the statistics of the data sample in the control re- 
gion. As such contaminants would raise the background 
estimate, we add a one-sided uncertainty in quadrature 
on the low side of the estimate, in amount equal to the 
statistical uncertainty in the estimate. 

The number of multijet events in the ep selection is es- 
timated by performing a fit of the C e distribution to tem- 
plates of true and fake electron taken from data. While 
the C e requirement must be relaxed for this fit, tight E/p 
requirements for the electrons and tight muon selection 
requirements are retained. The template distribution for 
true electrons is obtained from the ee pairs at the Z res- 
onance. The template distribution for fake electrons is 
obtained from like charge data events in a similar control 
region as is used in the ee and pp case. The method es- 
timates all contributions from processes containing fake 
electrons, including W + jet production, hence the cor- 
responding process is removed from the W MC sample. 

The shapes of the kinematic distributions for multijet 
backgrounds are modeled by events in like charge data 
with loose lepton quality but failing the tight criteria. 
The electron requirement is relaxed to C e > to obtain 
a sufficient number of events. The lepton quality require- 
ments are inverted for both leptons in ee and pp chan- 
nels. For the ep channel, events with a tight muon are 
also used to model the W + jet contribution included in 
the multijet rate prediction. The kinematic dependence 
due to the inversion of lepton qualities is corrected using 
the loose-to-tight efficiency for fake leptons obtained in 
the dijet sample and by the method described in the sec- 
tion for W boson background. For the pp selection, the 



multijet background mostly comes from the semilcptonic 
decays of heavy flavor quark decays; hence the correc- 
tion derived from the dijet sample is not used for the /x/x 
channel. 



MULTIVARIATE ANALYSIS 

A multivariate technique is employed to character- 
ize the Higgs boson signal and the backgrounds and to 
achieve maximum separation between them. A Boosted- 
Decision-Tree (BDT) algorithm [18] is used to construct 
a discriminant from kinematic variables taking into ac- 
count the variable correlations. The algorithm is robust 
against low number of events, which is particularly bene- 
ficial for the like charged dilcpton search where the data 
samples used to model the instrumental backgrounds are 
limited. 

The kinematic variables considered for the BDT inputs 
are: 

• Dilcpton kinematics: 

leading and trailing lepton transverse momenta 
[pfji, Vt\i invariant mass [M (£,£)}, angular sepa- 
ration [Ar](£,£), A(j)(£,£), AR(£,£)]- 

• Kinematics of all leptons in the event: 

lepton multiplicity [N ] , vector sum and scalar sum 
oi Pt of all leptons [p? 1 , ^Pt\'i 

• Kinematics of all jets in the event: 

jet multiplicity [AJ ot ], vector sum and scalar sum 
oi pt of all jets [Ht, Sp^ ]; 

• Kinematics of all the objects (leptons and jets) in 
the event: 

object multiplicity [iV ob J], vector and scalar sum of 
Pt of all objects [p T ° J , Hp T J ]; 

• Missing transverse energy: 

missing transverse energy \$t\ j component perpen- 
dicular to the object/muon which is closest to the 
$t axis in <f> {1£t spcc ^} to be insensitive to a pos- 
sible mismeasurement of the object momentum; 

• Dilcpton - T£t relation: 

transverse mass with minimum/maximum value 
[M T (^ T ) min / max ] (those calculated with re- 
spect to electron/muon for e/x channel only 
[Mx(e//x,^/r)]), minimum/maximum azimuthal 
angular separation [A0(^,$ T ) min / max ]. 

Between 11 and 17 variables are selected from this list for 
the training of BDTs based on the discrimination power 
of the variables for a given channel. Distributions of rep- 
resentative variables are shown in Figs. 1-3 for ee, /x/x 
and e/x channels after kinematic selection of like charge 
dilepton events. 



The BDT is trained for each Higgs boson mass consid- 
ered. The Higgs signal and background model samples 
described in the previous sections are separated into two 
orthogonal samples such that there is no overlap between 
events used for training and those used to derive the fi- 
nal result of the search. The training is carried out in 
two stages. The first stage uses one or two specific back- 
ground processes that are dominant after the kinematic 
selection; charge flip and multijet for ee and /x/x, multi- 
jet and W + jet for the e/x channel. The resulting in- 
strumental BDT discriminant, BDT^, is used to separate 
signal-like and background-like events (Fig. 4 a, c and 
c) . The second stage of training, referred to as "physics 
BDT" , uses those events that appear in the signal-like re- 
gion of BDTi , where the threshold is optimised for a sig- 
nal efficiency of approximately 90% for each mass point. 
The combined background samples are used for ee and 
/x/x training, which are mostly diboson and W boson or 
charge flip processes, and diboson only for e/x. All vari- 
ables from the list that have discrimination power and 
are well modeled by the prediction are selected for each 
stage of the training and for each of the three channels. 
The input variables used for the physics BDT training 
for ee and /x/x channels are: p^, pfp, Ar](£,£), M(£,£), 
p^ e , Y,p T , fl T , $t spoc , M T (£,# T ) m[n , A<^,$ T ) min / max . 
The input variables used for the physics BDT training 
for e/x channel are: p e ^, p^, Arj(£,£), A<j>(£,£), AR(£,£), 
M(£,£), p*, H T , 5V T et , p?° bj , Sp° bj , $ T , $ T spcc , 
M T (^,$ T ) min / max , A(/>(€,^T) min/max , BDT 4 . The final 
discriminant for each channel is computed as an effective 
product of the two discriminants after the selection based 
on BDT; (Fig. 4). 



CROSS SECTION LIMITS 

The final multivariate discriminants after all selection 
criteria (Fig. 4 b, d and f) show that the data are well 
described by the sum of the background predictions. In 
absence of an excess in the number of observed events 
over the SM backgrounds, upper limits on the production 
cross section have been determined. 

Uncertainties on the SM cross section for the pro- 
cesses modeled by MC simulation are 5% for the associ- 
ated Higgs boson production, 7% for diboson (WZ/ZZ) 
production, and 6% for W boson production. Experi- 
mental uncertainties assigned to the MC include a nor- 
malization uncertainty of 4.7%, which enfolds the lepton 
trigger and identification efficiencies and their kinematic 
dependences and the uncertainty of the Z boson pro- 
duction cross section used to study the normalization, 
and a 2% uncertainty on mismodcling of jets. An addi- 
tional uncertainty is assigned to the W boson background 
based on studies using the control samples; 20% uncer- 
tainty is used for each misidentified muon, 50%/14% for 
each misidentified electron originating from a photon/jet. 
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FIG. 1: (Color online) The distribution of (a) M(l,t), (b) AR(£,£), (c) JV j8t , (d) p£ obJ , (e) $ T 
channel comparing data and predicted backgrounds as well as the Higgs boson signal [Mh 
kinematic selection of like charge dilepton events. 



(f) A<P(e,^ T ) lnin , for the ee 
160 GeV) expectation after the 



The instrumental backgrounds, charge flip and multijet, 
which are estimated directly from data, have uncertain- 
ties between ll%-42%, arising mainly from the limited 
number of events in the control samples and the extrap- 
olation to the signal search region. The mismodeling of 
the kinematics of the instrumental backgrounds are rep- 
resented by a difference in the shape of the multivariate 
discriminants, and corresponds to up to 20% uncertainty 
on the final yield. The uncertainties described above are 
considered uncorrelated. 

The number of predicted and observed events after the 
kinematic selection of two like charged leptons and after 



the additional selection based on the multivariate dis- 
criminant BDTj are listed in Table I. The total uncer- 
tainties associated with each background and signal pro- 
cesses, excluding the shape uncertainty for the instru- 
mental backgrounds, are given. 

Cross section upper limits have been determined 
with the modified frequentist approach [19] with a log- 
likelihood ratio (LLR) test statistics. The systematic un- 
certainty of the signal and background predictions are 
represented by a Gaussian distributed fluctuation of the 
expected yield, where correlations across different chan- 
nels for any particular uncertainty are taken into account. 
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FIG. 2: (Color online) The distribution of (a) M(l,l), (b) AR(£,l), (c) N ict , (d) p^° bl , (e) ^ , (f) Ac/>(A^ T ) mil \ for the MM 
channel comparing data and predicted backgrounds as well as the Higgs boson signal (Mh = 160 GeV) expectation after the 
kinematic selection of like charge dilepton events. 



The details of the calculation are explained in Ref . [20] . 
Upper limits on the W/ZH cross section expressed as a 
ratio to the SM Higgs cross section and the correspond- 
ing LLR distributions are shown in Fig. 5. Ratio limits 
for each dilepton channel and for the combination are 
summarized in Table II. 



CONCLUSIONS 

A search for associated production of SM Higgs boson, 
pp — > W/ZH, has been performed with a final state with 



two like charged leptons, W/ZH -> ^ ± £' ± + X, in the 
ee, e/i and fifi channels. After all selection criteria, 17 
events in the ee channel, 24 events in the e/i channel, 
and 24 events in the \i[i channel have been observed in 
agreement with SM predictions. The observed (expected) 
upper limits on a(W/ZH) x B(W/ZH -> £ ± £' ± + X) 
for all three channels combined using a total integrated 
luminosity of 5.3 fb -1 collected by the DO detector in 
pp collisions at y/s = 1.96 TeV, expressed as a ratio to 
the SM Higgs cross section, are found to be 6.4 (7.3) 
for a Higgs boson mass of 165 GeV and 13.5 (19.8) for 
115 GeV. 
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FIG. 3: (Color online) The distribution of (a) M{£,£), (b) AR(£,£), (c) 7V jct , (d) p^ > , (e) 



(f) A<^, 



for the 



e/i channel comparing data and predicted backgrounds as well as the Higgs boson signal (Mh = 160 GeV) expectation after 
the kinematic selection of like charge dilepton events. Due to the absence of charge flip background in e/i channel, the sample 
after kinematic selection is dominated by multijet and W+jet backgrounds. 
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FIG. 4: (Color online) The distribution of (a, c, e) the multivariate discriminant against instrumental backgrounds, BDTi, and 
(b, d, f) the final discriminant which is an effective product of the BDT outputs from the first stage (instrumental) and the 
second stage (physics) for the (a, b) ee (c, d) fifi and (e, f) e/j, channels. Data and background predictions corresponding to 
an integrated luminosity of 5.3 fb -1 and the signal distributions are shown for an assumed Higgs boson mass of 160 GeV. The 
shaded region represents the total uncertainty on the background prediction. The final discriminant distributions are shown 
after requirements on BDT; > —0.85, —0.9, —0.2 for the ee, nfi, e/i channels, respectively, corresponding to approximately 90% 
signal efficiency points for a Higgs boson mass of 160 GeV. 
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TABLE I: The number of predicted and observed events for 5.3 fb" 1 of Run II integrated luminosity, after kinematic selection 
of like charged dileptons and after the final selection based on the multivariate discriminant against instrumental background, 
BDTi. This excludes the two control regions used for the estimation of the charge flip and multijet backgrounds. The 
uncertainties reflect the total systematic uncertainties. 
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TABLE II: The expected and observed production cross section limits in terms of ratios to the SM cross section for individual 
channels and for the combination. 
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FIG. 5: (Color online) (a) the expected and observed production cross section limits in terms of the ratio to the SM cross 
section as a function of the Higgs boson mass for an integrated luminosity of 5.3 fb -1 combining ee, fi/j, and e/j, channels, (b) 
the corresponding LLR distribution. 
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